Edge Occupation Numbers of a Half-Filled Landau Level 
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We demonstrate that a theory of the edge of a half-filled Landau level recently proposed by Lee 
and Wen predicts results for the edge occupation number similar to those of a variational trial wave 
function proposed previously by us. We treat Lee and Wen's edge action of a half-filled Landau level 
within the framework of bosonization theory, and show that the momentum occupation numbers 
are determined by a product of two Green's functions, one charged and one neutral. In the bulk 
region (k < 0) we find a linear occupation profile, nk <x A + Bk, while in the tail region (k > 0) it 
is exponentially decaying over the range k ~ A n , the momentum cutoff for neutral mode. We find 
a good fit with the numerical results for occupation numbers. 

PACS numbers: 73.20.Dx, 73.40.Hm 

Recently there has been a great deal of interests iriptwo-dimensional electron liquids in a strong magnetic field 
at half-filling, and several theories have been proposedErQ. Unlike insulating QH liquids this new state is metallic. 
Although all these theories describe a metallic state it is notrflear whether they predict identical edge properties. 
Field theoretical methods of Shytov et al.u and Lee and WenB seem to predict different dynamic properties of the 
edge. It is also not clear whether the edge properties of a variational trial wave function of this liquidQ differs from 
those predicted by field theories. 

After attaching two flux quanta to each electron to realize the Chern-Simons transformation at half-filling, one 
is apparently left with a Fermi sea without external magnetic field. One can alternatively regard it as N/2 of the 
composite fermions seeing an upward N/2 quanta of flux, while the other half sees magnetic flux in the opposite 
direction, in an iV-particle system. On average the net residual flux will be zero. As a consequence half of the 
composite fermions in the system occupy the "lowest Landau level" (as judged by the composite fermions), and the 
other half the anti-holomorphic lowest Landau level. This consideration naturally leadsrto a trial wave function at 
half- filling in particular, and wave functions at other hierarchy filling factors in generalliffl. 

Once a trial wave function for the half-filled Landau leveH is known the edge occupation number can be computed 
in the same manner as the calculation of the occupation number from Laughlin's wave function^. It turned out that 
the edge of a metallic liquid has a rather different property than that of incompressible QH states at filling factor 
v = 1/M (M odd integer). It is well known that for incompressible QH states at v = 1/M the electron occupation 
number i nsar the edge is oc (k e d — k) M ~ 1 0{k e d — k), where k e d is the wavenumber corresponding to the location 
of the edgeOO. Our result for the quantum Hall droplet at half filling showed that a tail exists in the occupation 
number beyond the radius of a droplet of uniform filling factor 1/2, i.e. ^ for k~k e d^. 0. This apparent increase 
in the size of the droplet is due to the displaced zeroes of wave function from the electronscE Thp_pccupation numbers 
were found to exhibit Friedel type oscillations, consistent with the metallic nature of the liquicfj. 

Recently Lee and Wen have proposed a twp-mode model for the edge! based on a theory of the electron liquid at 
v = 1/2 consisting of charge and flux liquidsu. In this theory the charge and the neutral modes of edge deformations 
are assumed to be well defined up to the momentum cutoffs A Cj „, beyond which modes become over-damped. Lee 
and Wen have calculated tunneling coefficient using this theory, and have found that it agrees with the experimental 
data of Grayson et alM 

In this paper we demonstrate that Lee and Wen's theory predicts the edge occupation number profile similar to 
what is found from the trial wave functionO. We treat their proposed action within the framework of bosonization 
theory and find that cx A + Bk for — A c <C k < and cx e~ fe / A " for k > (k is measured from k e d)- Our trial 
wave function also gives linearly increasing occupation numbers for small negative k. 

When a proper action is given for the edge at half-filing the electron Gxeen's function can be calclulated. The edge 
action of a half-filled Landau level derived from the dipole theory of Leea is 

* - h I £ (£-«£) ** « + / <*-»*•* - « 

It is important to note that the neutral mode is counter-propagating compared with the charged mode, as a result 
of the minus sign that comes in the dtr\d x r\ part of the action and the positive-definiteness of the Hamiltonian. The 
last term in Eq. (Q) is the coupling of the edge mode with the transverse bulk current (jt-n) or the tunneling (d^j^) 
current. We parameterize the sample boundary with x, using n to indicate the outward normal to the boundary. We 
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quantize the edge action in the absence of the coupling to current sources. The canonical commutation relations for 
the fields <fi and rj are 

&-,<t>{v)\ = m5(x - y); [^,r]{y)} = -2mS(x - y). (2) 



By comparison with the chiral Luttinger liquid theory at other filling factoraia'ta, one may associate 2d<p/dx with the 
charge displacement at the edge of the liquid. Since the one-dimensional density is (v/2ir) (displacement), electron 
creation/annihilation operators at the boundary will need to satisfy the following algebra: 

[—(displacement), (electron op.)] = ±S(x — y)(electron op.) (3) 

One can easily show that e ±2i * satisfies the desired algebra, but when one interchanges two such operators at x and 
x', one finds that they obey Bose, rather than Fermi, statistics. To fix the statistics without altering the relation in 
Eq. (|J), we use the following annihilation/creation operators 

In Lee and Wen's paper, a specific combination 2<f> — rj appears in the electron Green's function, while our formulation 
seem to favor both signs. One presumably needs the coupling term in Eq. ([!]) in order to discriminate the two signs. 
It will be clear that all the computations below will be insensitive to the choice of sign. We introduce the mode 
expansion of the boson fields: 

<j>{x) = fo+fox+T ^=[ale- mx -a n e mx ]- [f , M = i/2, 

n>0 ^ 2U 



T)(x) = m+9ox+ V-^[6„e lnx -ble mx ]- [30,770] = -1. (5) 

n>0 v 



The ground state is defined by the action of the annihilation operators, a n \gs) = b n \gs) — 0. The time dependence is 
given by setting x — > x — v c t for the charged mode, and x — ► x + v n t for the neutral mode. 

Since this is a bosonization theory, the phonon Green's functions determine the underlying electron Green's function 
as well. Define z = x — v c t, and w = x + v n t. We have 

roc ii 

G^z) = 4(0(z)0(O)) - 4(0(0)0(0)) = 2 / — [e ikz - l] e - fe / A = = -21n[l - iA c z], 

Jo k 

G v (w) = (r)Mv(0)) - <»?(0)»j(0)> = - ln[l + iA n w]. (6) 
Electron Green's functions are obtained as 

$(*M0)) = ( e a*{«) e -M«o))( e -*?W e *i(o)) = exp[G (z) + G n (w)], 

(^(0)^(z)) = (e-^ )^')) (e^e^W) = exp[G (-z) + G v (-w)]. (7) 
It also follows from the definition in Eq. (0) that (ip(z)'ip(O)} = (^j(z)ijj(0)) . Putting our results together, 

(il>(z)${0)) = (Hz)iP(O)) = 1/[1 - «A c z] 2 [l + iA n w], 

(i/)(0)1>(z)) = ($(p)il>(z)) = 1/[1 + iA c z] 2 [l - iA n w}. (8) 
The time-ordered Green's function G e (x,t) — 9(t)(ijj(z)ip(0)) — 6(—t)(ip(0)iJ;(z)) is given by 

e[ ' ' [1 - iA c z] 2 [l + iA n w] [l+iA c z] 2 [l-iA n w]' [ 1 

One can consider another kind of Green's function, Gb(x, t) = exp[8(t){G,f ) (z) + G n (w)} + #(— rj){G^(— z) + G n (—w)}}: 

r( t\= e[t) 1 °k*l nm 

h[ ' ' [1 - iA c z] 2 [l + zA„ w ] [1 + iA c zY[l - iA n w] ' 1 ' 

This function is clearly bosonic in character, Gb(— x, — t) — +Gb(x,t). On the other hand, the electron Green's 
function satisfies G e (—x,—t) = —G e (x,t). Notice that one has the approximate symmetry G e {x,—i) ~ —G e (x,t) 
provided A c x, A n v n t <C 1, because we have, in such a limit, 
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We will assume that the momentum cutoffs satisfy A„ « A c . For the bosonic Green's function Gb(x, t), one observes 
a similar symmetry in the slightly different limit A n v n t Cl< A n x, and x -C v c t; 

G b (x, t) « ^ I — - — ) = « -G t (g, -t). (12) 

The above expression embodies the symmetry that was extensively discussed and emphasized by Lee and Weni. It 
says that the desired statistics property shows up at distances short compared with v c t, but larger than 1/A„. 

We now consider the static property of the edge, in particular its momentum occupation. The momentum distri- 
bution is simply a Fourier transform of the equal-time Green's function G e (x,Q + ) = Gb(x,0 + ): 



n k oc 



oo ikx J T A A 

.„ [1 - tA c x] 2 [l + iA n x] y J A n + A c w v ' 



Due to the presence of counter-propagating neutral mode, we have an extra occupation on the k > side that decays 
exponentially over the distance k ~ A„. We believe the discontinuity at k = is not physical and is an artifact of 
taking the integration over infinite range. We can instead work with the Green's function for finite circumference L: 



G^z) = -2 In 



l-exp[-27r(l-z;zA c )/A c L] 
l-exp[-27r/A c L] 



= _ lll( l-exp[-2 7 r(l+mA ra )/A w L]- 
1 — exp|— 2tt/ A n L\ ' 



The corresponding occupation number n m is given by 



n m oc 



exp[27rima;] dx 



[1 — exp(— 27r(l — in c x)/K c )] 2 [l— exp(— 2Tr(l + iK n x) / n n )] 



(15) 



We have introduced k c .„ = A C .„L and measured x in units of L. Expanding the denominators in powers series, one 
obtains a simple expression for n m : 



oo oo 

II n, X 

s=0 t=0 



£ E S ts,rn(s + l) e - 2 -/«c- 27 r t / Kn _ (lg) 

Normalizing by the value at m = 0, 

_f e -W«n i m>0 
n m /?i - | e 2. m / Kc{1 „ m(1 „ e -2x( Kr 1 + «- 1 ) )}] m < o l 1 ^ 

It is exponentially decaying in the m > region (tail) and linear, of the form A + Bm, for — k c <C m < 0. The 
coefficients A and B depend on the exact value of the cutoff parameters. To connect with the thermodynamic limit, 
one can rewrite the occupation numbers in terms of k. We assume that the momentum cut-offs A n and A c remain 
finite in the L— >oo limit. 

Uk/no ^ { eX ( i_ fe (A-i+A-i)), I < 0. (18) 

Figure 1 shows occupation numbers for N = 11 electrons at half- filling obtained from finite-size exact 
diagonalizationcl. Our trial wave function indicates that the tail region, x > with x — kR e d, does not disap- 
pear in the thermodynamic limit. Here R ec i — 2\/~N is the radius of a uniform, compact half-filled droplet. The 
trial wave function was constructed by occupying about half of the composite fermions in the lowest Landau level, 
while each of the other half singly occupies the higher Landau levels after the projection. In— the Laughlin state all 
of the composite fermions occupy the lowest Landau level to form an incompressible liquidu-3. Such differences in 
the occupation behavior of composite fermions underlie the observed difference of electron occupation for compress- 
ible/incompressible states. From Eqs. ( |l7| ) and Jllj ) we expect n k = A + Bk for small negative k, and we find that 
this expression (open squares) fits the numerical results well for suitable choice of A and B. 
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In this paper we have demonstrated that the edge occupation numbers calculated from our trial wave function are 
similar to those of Lee and Wen's theory. It would be interesting to see if similar excitation properties can be found. 
The work of SREY has been supported by the Ministry of Education under Grant No. BSRI-96-2444. 
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FIG. 1. Occupation numbers for N = 11 
fitted with a linear function (open squares). 
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computed numerically. Occupation numbers for small negative values of x are 
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